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Abstract 

We present an extraction of azimuthal correlations between two pairs of charged pions detected in 
opposite jets from electron-positron annihilation. These correlations may arise from the dependence 
of the di-pion fragmentation on the polarization of the parent quark in the process e'’“e“ —)• qq. 
Due to the correlation of the quark polarizations, the cross-section of di-pion pair production, in 
which the pion pairs are detected in opposite jets in a dijet event, exhibits a modulation in the 
azimuthal angles of the planes containing the hadron pairs with respect to the production plane. 
The measurement of this modulation allows access to combinations of fragmentation functions that 
are sensitive to the quark’s transverse polarization and helicity. Within our uncertainties we do 
not observe a significant signal from the previously unmeasured helicity dependent fragmentation 
function Gi . This measurement uses a dataset of 938 fb“ ^ collected by the Belle experiment at or 
near ^/s ^ 10.58 GeV. 

PACS numbers: 13.88.-l-e,13.66.-a,14.65.-q,14.20.-c 
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Fragmentation functions are fundamental QCD objects describing the hadronization of 
quarks into final-state hadrons. They cannot be computed using lattice QCD techniques 
due to the non-inclusive hnal-state and therefore have to be measured. In addition to 
studying properties of QCD, fragmentation functions are important tools to extract the 
partonic spin structure of the nucleon from semi-inclusive deep-inelastic scattering (SIDIS) 
and proton-proton (p-l-p) collision data [H [2] . In recent years the role of polarization depen¬ 
dent fragmentation functions as quark polarimeters has become more important to access 
polarization dependent quark distributions in SIDIS and p-|-p [3]. 

Di-hadron fragmentation, e.g., q —)■ {h^h~) + X, describes the hadronization of a quark 
into a hnal-state containing two distinguishable hadrons. Compared to single-hadron frag¬ 
mentation, it offers access to transverse polarization of the fragmenting quark without in¬ 
volving either hnal-state polarimetry or transverse-momentum dependence of the fragmen¬ 
tation process. This is achieved through the additional degree of freedom provided by the 
relative momentum of the hadron pair. In case of the chiral-odd di-hadron fragmentation 
function the orientation of the plane spanned by this relative momentum and the total 
di-hadron momentum can be correlated with the transverse polarization of the fragment¬ 
ing quark. Such correlation is reminiscent of the explicitly transverse-momentum dependent 
Collins fragmentation function where the transverse momentum of the outgoing hadron 
(or di-hadron) provides the plane correlated with the quark polarization. However, such cor¬ 
relation and thus does not survive simple integration over transverse momentum, which 
complicates the theoretical description. In contrast, the di-hadron plane does survive such 
integration, permitting application of the ordinary collinear framework. 

The chiral-odd has already been subject to measurements in e''“e“ annihilation into 
di-hadrons |1] as well as in SIDIS and p-|-p. Its chiral-even partner, the helicity-dependent 
di-hadron fragmentation function Gi, has yet to be measured. Together with the chiral-odd 
fragmentation functions and the quark-polarization independent Di, it completes the set 
of leading-twist di-hadron fragmentation functions. Unlike Hf, which has a counterpart in 
the single-hadron fragmentation with the Collins fragmentation function, for Gj- no such 
correspondence exist [3]. This makes it an elusive tool to probe quark helicity (e.g., in the 
nucleon) through fragmentation. In addition, it has been conjectured that the measurement 
of G;j*“ in e+e“ —)■ {hih 2 )(hih 2 ) + X, i-e., production of two di-hadrons from both the quark 
and the anti-quark in e’''e“ annihilation, is related to jet-handedness correlations, which 
might receive contributions from GP violating effects [5]. 

The process of di-hadron production from both the quark and the anti-quark in 
annihilation is the cleanest way to access polarization-dependent di-hadron fragmentation. 
While in SIDIS and p-|-p the fragmentation functions are folded with partially unknown 
parton distributions, such “contamination” is missing in e+e“ annihilation. Furthermore, 
the spin-correlation of the quark anti-quark pair in e’''e“ annihilation allows the study of 
polarization dependent fragmentation functions even in the case of an unpolarized initial 
state. For the study of —>■ {hih 2 ){hih 2 ) + X has another advantage: formally G^ 

is a transverse-momentum dependent fragmentation function and as such does not survive 
simple integration over the transverse momentum generated in the fragmentation process. 
However, the correlation between the two hadron pairs in this specihc process introduces 
a particular weighting factor in the integration, such that the cross section itself survives 
integration over the observed transverse momentum. The fragmentation functions will thus 
enter as specihc transverse-momentum moments in the cross section expressions. However, 
we shall adopt the simplihed notation of G^{z,M‘^) from Ref. [3], where 2 ; is the fractional 
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momentum carried by the hadron pair and M its invariant mass. 

This paper reports on the first measurement of the G;j*--induced correlation between two 
di-pions in the process e’''e“ —?• (7r’''7r“)(7r’''7r“) + X, where the two di-pions are from two 
separate and opposing jets. 

The Belle experiment [B] at the KEKB storage ring recorded about 1 ab“^ of 
annihilation data. The data was taken mainly at the T(4S') resonance at ^/s = 10.58 GeV, 
but also at T(IS') to T(5S') resonances and at ^/s = 10.52 GeV. The Belle instrumentation 
used in this analysis includes a central drift chamber and a silicon vertex detector, which pro¬ 
vides precision tracking for tracks in 17.0° < 0Lab < 150.0°, and electromagnetic calorimeters 
(EGL) covering the same region. Particle identihcation is performed using information on 
dE/dx in the GDG, a time-of-flight system in the barrel, aerogel Gherenkov counters in the 
barrel and the forward endcap, as well as a muon and identihcation system outside the 
superconducting solenoid, which provides a 1.5 T magnetic held. 

Using the event shape variable thrust, a sample of light and charm quark {u,d,s,c) 
fragmentation can be selected HlIH]. The thrust is dehned as 


T 


max 

|n|=l 


Ei 1^ • Pi 

Y.i \Pi\ 


( 1 ) 


Here the pj are all particles in the event. We use a thrust cut of T > 0.8 which reduces the 
contribution of B mesons to less than 1%. Since the polarization-dependence of the fragmen¬ 
tation would average out in unpolarized annihilation into one di-pion, a measurement 
has to use the correlation between the polarizations of the produced quark-antiquark pairs. 
This then leads to a correlation of hadron-pairs detected in opposite hemispheres [3]. The 
correlation is expressed as a modulation of the di-hadron pair cross-section in the angles 
and They are defined in Ref. |3] as the azimuthal angles of the planes containing the 
hadron pairs with respect to the event plane spanned by the beam axis and the momentum 
sum vector Ph = hi + h 2 of one hadron pair in the Genter-of-Mass System(GMS), where hi, 
^2 are the momenta of the hadrons forming the pair. It should be noted, that the reference 
frame defined here, is asymmetric with respect to the hadron pairs, both angles and 

are measured around Ph- With this dehnition, the di-pion fragmentation function Gj; 
can be accessed in the moment 


(cos(2(«hR; - cx (2) 

Here q, q are the quark and antiquark flavors, Cq the corresponding electric charges and 
M the fractional energy and invariant mass of the hadron pair produced in the antiquark 
fragmentation. 

Instead of using the asymmetric frame defined above, this analysis uses a symmetric 
frame in which the angles are computed using the jet axes of di-jet events. The jet axes 
are reconstructed using the anti-fc^ jet algorithm [S]. They serve here as proxies of the 
directions of the fragmenting (anti)quarks. A similar approach had been employed in the 
previous measurements of and the Gollins fragmentation function |H [8] in a symmetric 
frame using the thrust axis. The coordinate system is depicted in Fig. [1} 

The angles <I)r. can be expressed as: 

^Ri = sgn (ji - ((ex ji) X {ji X Ri)')') arccos [ . (3) 

^ ^ Vlexjil \jiXRi\J 
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FIG. 1: The coordinate system used for this measurement. For convenience, the jets are depicted 
to he in the event plane which in turn is spanned by the beam axis and the jets. 


Here e is the e -beam direction, ji the jet axes and Ri = h\ — with h\ the momentum 
of the positive pion and momentum of the negative pion in jet i. All quantities 

are in the CMS. Note that one of the jet-axis directions is flipped such that in an ideal 
back-to-back event, both jet axes point in the same direction. Compared to the asymmetric 
frame, the theory for the symmetric frame has not been worked out completely, though we 
believe it to be just a modihcation of Eq. ([^. Compared to |1] and jH], where the thrust 
axis of the event was used as a proxy for the quark direction, the jet method will be more 
robust against hnal-state gluon radiation, which is usually described by Sudakov factors or 
evolution equations. 

We use jets reconstructed with the anti-/cr algorithm as implemented in the FastJet 
3.0.5 package [TOl n] with a jet radius of R = 1.0 and a minimum energy of 3.75 GeV 
for each jet. This cut on the jet energy also effectively removes e’''e“ —>■ rr events and 
leaves only dijet events. We restrict the allowed range of the polar angle of the jets in 
the CMS to 1.38 rad < ^cms < 1-75 rad to minimize systematic effects at the edges of 
the barrel acceptance. Jets are reconstructed from all charged tracks that originate from a 
region around the reconstructed interaction point (IP) as well as photons with energies larger 
than 50 MeV in the barrel, 75 MeV in the forward endcap and 100 MeV in the backward 
endcap. The IP is required to fulhll \dr\ < 2 cm and \dz\ <4 cm, where dr and dz are the 
distance of closest approach to the interaction point in the plane perpendicular to the beam 
direction and along the direction of the beams. We select charged pion pairs by using all 
charged tracks, which have not been identihed as kaons, protons, electrons, or muons [12] 
and requiring a minimum fractional energy z of the individual tracks of 0.1. 


We perform a two-dimensional fit fo fh® normalized di-pion pair yields 


(JV> 






FIG. 2: Results for binned in M and 2 ;. The black error bars are statistical and 

the green bands show the systematic uncertainty. 


using the truncated Fourier expansion 


1 + cos + ^ R ,) + . 


(4) 


The amplitude '*’« 2 )) corresponds then to the quantities in Eq. (2), whereas the 

amplitude corresponds to the previously published measurement of |3]. 

Results binned in M and 2 ; of one hadron pair are shown in Fig. [^integrated over 2 ; and M 
of the other pair. 

The purity of the sample is above 90% for di-pion pairs, with the remaining 10% com¬ 
ing from an admixture of kaons. The leading systematic effects are the contributions from 
detector smearing and hducial cut edge effects to the measured asymmetries. These are 
estimated from simulations using Pythia [13] and EvtGen na for various physics processes 
not including the polarization dependent ones of Eq. (j^, and GEANT3 [T5| for the detector 
effects. In addition to the light {u, d, s) quark fragmentation, our results also contain signih- 
cant contributions from charm quark production. The charm contribution is 18% on average 
and approximately constant over the invariant mass bins. In z, it drops monotonically from 
21% to 5%. 

Within our statistical and systematic uncertainties, we do not observe a deviation from 
zero for Statistical uncertainties are symmetric and about 10“^ in mag¬ 

nitude, whereas our systematic uncertainty is asymmetric and varies between 2 ■ 10“^ at 
low M and low ^ and 5 ■ 10“^ at high M and z. It has been conjectured that the Gi 
fragmentation function receives non-zero contributions mainly from p-p wave interference in 
the partial-wave expansion of We note that an enhancement of such contribution by 
restricting 6*decay, fhe decay angle of the positive hadron in the hadron-pair GMS used in 
such expansion, to values for which cos(6*decay) is positive does not result in a signihcantly 
non-zero amplitude either. 

In summary, we show hrst results for azimuthal modulations in the cross-section of di- 
pion pairs in di-jet production from electron-positron annihilation. The amplitude of the 
cos(2(<I)i^j — $^ 2 )) modulation, which is sensitive to the helicity dependent fragmentation 
function G^, is consistent with zero within our statistical and systematic uncertainties. 
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